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Chiral and Critical Behavior in Strong Coupling QCD 
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We use a cluster algorithm to study the critical behavior of strongly coupled lattice QCD in the chiral limit. 
We show that the finite temperature chiral phase transition belongs to the 0(2) universality class as expected. 
When we compute the finite size effects of the chiral susceptibility in the low temperature phase close to the 
transition, we find clear evidence for chiral singularities predicted by chiral perturbation theory (ChPT). On the 
other hand it is difficult to reconcile the quark mass dependence of various quantities near the chiral limit with 
ChPT. 



1. INTRODUCTION 

One of the important challenges in lattice QCD 
is to compute quantities that are dominated by 
the physics of light quarks. Although there has 
been substantial effort in extracting such quanti- 
ties from the lattice results by matching the data 
with ChPT, it is unclear if the current range of 
masses used in the calculations are in the range 
where ChPT is valid 1. Recently, efforts have 
been directed in two directions: a) Improving 
actions that reduce lattice artifacts, b) Improv- 
ing ChPT that take these artifacts into account. 
However, we think that it is also equally impor- 
tant to find improved algorithms to approach the 
chiral limit. 

Over the last decade a new class of cluster al- 
gorithms have emerged for solving a variety of 
lattice field theories 2 . These algorithms help in 
beating critical slowing down very efficiently. Re- 
cently, it was shown that the strong coupling limit 
of lattice gauge theories with staggered fermions 
can be solved using these new algorithms |3] . For 
the first time, this allows us to study the physics 
of massless quarks on large lattices from first prin- 
ciples. 

In this article we present results from our study 
of the chiral and critical behavior near the chiral 
phase transition in strongly coupled lattice QCD 



with staggered fermions. We use £7(3) gauge 
fields instead of SU (3) in order to avoid technical 
complications in the algorithm. Further details 
of our study can be found in 4 . 

2. MODEL AND OBSERVABLES 

The model we study can be specified by the 
Euclidean action, 



-m^2 jjxi>x 



(1) 
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where m is the staggered quark mass, U x<tl are 
U (3) link variables and ip, ip are Grassmann vari- 
ables representing the staggered quark fields. We 
choose the staggered fermion phase factors r) x ,fj, 
to have the property that rfc ^ = l,fJ, = 1,2,3 
(spatial directions) and rft. 4 = T (temporal di- 
rection), where the real parameter T acts like a 
temperature. By working on asymmetric lattices 
with L t << L and allowing T to vary continu- 
ously, one can study the finite temperature phase 
transition in strong coupling QCD 0. 

The model described by the partition function 
Z(T,m), constructed from the above action in 
the usual way, is known to have an exact 0(2) 
chiral symmetry when m — 0. This symmetry 
is broken at low temperatures but gets restored 
at high temperatures due to a finite temperature 
chiral phase transition. In order to study the chi- 
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ral physics near this transition we focus on the 
chiral condensate, 



the chiral susceptibility, 

i i a 2 
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-Z{T,m), 



L 3 Z dm 2 
and the helicity modulus, 
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(4) 



N/8), with a x = 1 on 
even sites and a x = — 1 on odd sites. When 
m = the current J X4Ji is the conserved current 
associated with the 0(2) chiral symmetry. Fur- 
ther, as discussed in 6 , it can be shown that 
F 2 = limi^oo Y m (m ~ 0), where F is related to 
the pion decay constant. We also measure the 
Goldstone pion mass M„. 

3. UNIVERSAL PREDICTIONS 

The predictions of ChPT for 0(N) models have 
been discussed in [Sj. In particular the finite size 
scaling formula for \ at m = is given by 



X = -S 2 L 3 
N 



, (5) 



where N = 2 in our case, (3% — 0.226..., E = 
lim m ^o hnii^oo^), and a is a constant depen- 
dent on other low energy constants. The quark 
mass dependence of (</>), x an d Y m are given by 



<</>) 



1 + u\\pm + a2r 



Y m — F 2 [1 + a 3 \fm + a±m 
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M 2 = — -r-[l + a 5 ^/m + a 6 m + 



(G) 
(7) 

(8) 



For N = 2 one further finds that = 0. One 
of the consequences of chiral symmetry can be 
described by the relation 



(4)i 



M 2 F 2 

7T m 



1 + C(m), 



(9) 



If the chiral phase transition is second order 
then close to the critical temperature T c we have 

E(T) = A(T C - Tf, T < T c , (10) 
lim ((f)) = Bm 1/S , T = T C , (11) 

L — >oo 

The 0(2) universality predicts predicts (3 = 
0.3485(2), 8 = 4.780(2) 

4. RESULTS 

We have done extensive calculations on various 
lattice sizes in the range 8 < L < 192 with L t = 4. 
In order to understand the critical behavior we 
have measured \ f° r m = at different values 
of T between 7.3 and 7.5. On the other hand 
we focused on a single temperature in the broken 
phase (T = 7.42) and computed (</>), Y m and M v 
for masses in the range < m < 0.01. In figure 
n we plot our results for % as a function of L 
at m = and T = 7.42. The data fits well to 
the ChPT prediction (eq.0) with E = 1.079(2), 
F = 0.181(4) and a = 114(4) with a x 2 / d -0-f of 
0.73. The value of F obtained from this fit is in 
excellent agreement with F = 0.181(1) obtained 
through a direct evaluation of Y m at m — at 
large volumes. This can be seen from the plot 
shown in the inset of figure ^ 
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Figure 1. Plot of x vs - L an d vs. L (inset) at 
T = 7.42 and m = 0. 



which is the Gcllmann-Oaks-Renner relation. 
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Figure 2. Plots of (</>) (upper-left) Y m (upper- 
right), Mtt (lower-left), extrapolated to L = oo, 
as a function of the quark mass. Check of the 
Gclhnann-Oaks-Renner relation (bottom-right). 



Given this excellent agreement with ChPT at 
7tl = we have also looked for the quark mass 
dependence of ((f)), Y m and M T . These quanti- 
ties were computed for various volumes until their 
values did not change for two different volumes. 
This "infinite" volume data is plotted on figure^ 
When the data is fit to the prediction from ChPT 
(eqs. ©,0 and ©), with S and F fixed to 1.079 
and 0.181 (see above), we find that ai — 17.9(2), 
a 2 = -64(3), a 3 = 20(1), a 4 = 301(36), a 5 = 
-8.14(10) and a 6 = 38(2). The x 2 /d.o.f of all the 
fits are close to one. Although these fits appear 
to be good, the results are not consistent with 
ChPT. In particular ChPT predicts a 3 = 0. 
In spite of this disagreement, the bottom right 
plot of figure |21 shows that our results satisfy the 
the Gelhnann-Oaks-Renner relation (eq.(©). 

The only plausible explanation we can imagine 
for the disagreement, is that most of the masses 
used in the fit are perhaps too heavy for ChPT 
to be valid. This is because the cutoff in ChPT 
is proportional to F 2 (we are effectively in three 
dimensions), which is quite small in our case. In- 
deed the coefficients on are also uncomfortably 
large for the same reason. We are currently study- 
ing another temperatures deeper inside the bro- 
ken phase where F 2 is larger. 
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Figure 3. Plots of S vs. T and (<fi) vs. m at 
T = T c (inset). 



Finally a careful analysis of the finite size ef- 
fects of x & t various temperatures reveals that 
T c = 7.47739(3). We also find that the data for £ 
and ((f)) fits well to the form predicted by eqs. I jlOjl 
and G3J. We find (3 = 0.348(2), A = 2.92(2) 
with a x 2 /d.o.f of 0.53 for the fit to eq. ijTUj) and 
6 = 4.97(10), B = 5.9(1) with a x 2 /d.o.f of 0.34 
for the fit to eq. iJUJ. The results and the fits are 
shown in figure 
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